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Abstract. We consider the effect of Coulomb interactions on the average density of states (DOS) of dis- 
ordered low-dimensional metals for temperatures T and frequencies ui smaller than the inverse elastic 
life-time I/tq. Using the fact that long-range Coulomb interactions in two dimensions {2d) generate In^- 
singularities in the DOS vii^i) but only In-singularities in the conductivity o-(tj), we can re-sum the most 
singular contributions to the average DOS via a simple gauge-transformation. If limi^^o o"(tj) > 0, then 
a metallic Coulomb gap v{lo) oc \(j->\/e* appears in the DOS at T = for frequencies below a certain 
crossover frequency il2 which depends on the value of the DC conductivity (t(0). Here, — e is the charge 
of the electron. Naively adopting the same procedure to calculate the DOS in quasi Id metals, we find 
v{uj) oc (|tj|/J7i)^''^ exp( — I7i/|aj|) at T = 0, where J7i is some interaction-dependent frequency scale. How- 
ever, we argue that in quasi Id the above gauge-transformation method is on less firm grounds than in 2d. 
We also discuss the behavior of the DOS at finite temperatures and give numerical results for the expected 
tunneling conductance that can be compared with experiments. 

PACS. 71.10.Pm Fermions in reduced dimensions - 71.23.-k Electronic structure of disordered solids 
- 71.30.-|-h Metal-insulator transitions and other electronic transitions - 72.15.Rn Localization effects 
(Anderson or weak localization) 



1 Introduction 

In the early eighties, Altshuler and Aronov [Q perturba- 
tively studied the effect of electron-electron interactions 
on the density of states (DOS) of low-dimensional weakly 
disordered interacting electronic systems. For tempera- 
tures T and frequencies lo smaller than the inverse elas- 
tic life-time I/tq, they found that in reduced dimensions 
the interplay between disorder and electron-electron inter- 
actions gives rise to singular corrections to the averaged 
DOS. In two dimensions (2d), the long-range Coulomb in- 
teraction gives rise to a In^-correction to the average DOS 

El, 



Fermi energy (per spin projection) of electrons with effec- 
tive mass m and charge — e, Dq = vpi/2 is the diffusion 
coefficient in 2d, and ^ = wfTq is the elastic mean free 
path. We use units such that h ^ ks = ^- Note that for a 
good metal at high densities the Thomas-Fermi screening 
length is short compared with the elastic mean free path 
(k£ S> 1) so that Ti <C tq. In quasi Id metallic wires (which 
consist of many transverse channels but permit diffusive 
motion only in one direction) the leading correction to the 
average DOS is 



v(uj) 



where 



l-^ln(|c.|ro)ln(|c.|ri)' 



2e' 



■To 



{27T)^iyoDo (27r)2(To nkpi 



2d , (1) 



(2) 



7r|w| 



quasi Id , 



(4) 



where vq — 1/ {t^vf) is the DOS per spin in Id, and the fre- 
quency scale /2i depends on the effective electron-electron 
interaction constant /o. 



is a dimensionless measure for the resistance of the system 
at frequency scale w « Tg~^ (where (Tq is the Drude con- 
ductivity) , and the interaction-dependent time ri is given 

by 

Here, k — 2me^ is the Thomas-Fermi screening wave- 
vector in two dimensions, i/q — m/2TT is the DOS at the 



ill 



32ttDq 



(5) 



Here, Dq = vpi the (bare) diffusion coefficient in quasi 
Id. 

Obviously, the correction terms in Eqs. (|l|) and (^) 
diverge for w — > 0, so that at low frequencies these per- 
turbative expressions cease to be valid. What is the true 
low-frequency behavior of the DOS of disordered metals in 
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reduced dimensions? The answer to this question is rele- 
vant for a number of recent tunneling experiments 
where a strong suppression of the tunneling conductance 
G{V) as a function of the applied voltage has been ob- 
served (zero bias anomaly) . The tunneling conductance is 
related to the DOS via G{V) oc v{lo — eV), so that the ex- 
perimentally observed zero bias anomaly in the tunneling 
conductance reflects the strong suppression of the average 
DOS at the Fermi energy. 

In 2d, the low-temperature behavior of the DOS of 
a strongly correlated disordered metal has recently been 
measured by Bielejec et al. While at higher temper- 
atures they found a logarithmic correction, at the lowest 
temperatures they found a stronger-than linear suppres- 
sion of the DOS, which has been interpreted in terms of 
a hard correlation gap. The knowledge of the low-energy 
behavior of the DOS of a 2d disordered metal with a fi- 
nite conductivity might also be important to gain a better 
understanding of the physical mechanism that is respon- 
sible for the apparent metal-insulator transition in doped 
semiconductor devices [||,|6|. An intensely studied quasi 
Id system where under certain conditions the electrons 
propagate diffusively in only one direction are multi-wall 
Carbon nanotubes P,^. 

Let us briefly review previous calculations of the zero 
bias anomaly. In 2d, the first attempt to determine the 
true low-frequency asymptotics of v{uj) was apparently 
due to Finkelstein |^] who found that 1^(0-") oc for 
w — > 0. However, in the derivation of this result he as- 
sumed that the conductivity (t{uj) diverges logarithmically 
for w ^ 0. The behavior of the DOS of 2d disordered 
electrons with a finite conductivity was not calculated by 
Finkelstein. Later this problem was reconsidered by Be- 
litz and Kirkpatrick j^, who found that for frequencies 
exceeding the crossover frequency 



= Tg ^ exp[-l/ro] , 



(6) 



the perturbative expression given in Eq. (|I|) can actually 
be exponentiated, so that 



v{u}) ~ vq exp 



ln(|tj|ro)ln(|w|Ti) 



(7) 



This expression has been re-derived in different ways by 
several authors Note, however, that Eq. (0) is 

not valid for frequencies smaller than the crossover fre- 
quency ^2 defined in Eq. (|) |,|l^,|ri). A simple inter- 
polation formula for the DOS, which yields a physically 
sensible result even for u; ^ 0, has been proposed by one 
of us in Ref . [0 . This formula is based on a re-summation 
of the leading In^ -singularities to all orders in perturbation 
theory, consistently neglecting sub-leading terms that in- 
volve only logarithmic corrections. In this approximation, 
one obtains 



v{uj) « vo- 



' sin[^ 
dt exp 



\n{t/Ta)Ht/Ti) 



(8) 

For \bj\ ^ = Tq^ exp[— 1/ro], this expression reduces to 
Eq. (0). Note that Eq. (||) amounts to an exponentiation 



of the perturbative result in the time domain, whereas 
in Eq. (7|) the perturbation series is exponentiated in fre- 
quency space. In the limit cj — > 0, the leading term in the 
asymptotic expansion of Eq. (H) is 



.1/2 



1/2 



H 
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(9) 



Noting that ^toTi = ro/(j^oe'*), Eq. (||) can also be written 
as 

Hu^) - , (10) 
where the numerical constant C is given by 

C7 = 4(ro/7r)i/2exp[l/ro] . (11) 



Surprisingly, Eq. ( |10[) resembles the well-known classical 
Efros-Shklovskii Coulomb gap of two-dimensional elec- 
trons in the localized regime, where the DC conductivity 
(t(0) vanishes Note, however, that in the derivation of 
Eq. ( |l0| ) we have assumed that the DC conductivity cr(0) 
remains finite, i.e. the electrons are assumed to be delocal- 
ized. The dilference between the localized and the delocal- 
ized regime manifests itself in the dimensionless prefactor 
C: whereas in the case of the classical Coulomb gap of lo- 
calized electrons the constant C is a number of the order 
of unity that depends on the geometry of the underly- 
ing lattice, for the quantum Coulomb gap discussed here, 
C depends on the dimensionless conductivity of the sys- 
tem. The existence of the Coulomb gap in the delocalized 
regime of a disordered interacting 2d metallic system was 
also found numerically by Efros and Pikus (l3|. More re- 
cently, an intermediate delocalized phase in small clusters 
of disordered interacting electrons has been found numer- 
ically in Ref. (w). 

In deriving Eq. (pd]), we have implicitly assumed that 
the DC conductivity (t(0) does not deviate significantly 
from the conductivity <j(uj) at frequency w ~ I/tq. We 
have argued elsewhere that a finite renormalization of 
the conductivity can be approximately taken into account 
be replacing rp r» in Eq. ( |Tl| ) where r» = e^/(27r^cr(0)) 
corresponds to the true DC conductivity of the system. 
The connection between low-frequency behavior of the 
conductivity and the DOS has also been emphasized by 
Nazarov , and by Levitov and Shytov |^ . 

The zero bias anomaly in Id has so far received much 
less attention than the corresponding anomaly in 2d. Re- 
cent experiments on metallic Carbon nanotubes have mo- 
tivated Mishchenko et al. to study the fate of the 
perturbative Altshuler-Aronov correction in Id at low fre- 
quencies. They found that long-range Coulomb interac- 
tions in a quasi Id metal lead for — > to an exponen- 
tially small DOS, 



^{uj) oc exp 



^^InfA 



(12) 



where eo and ei are some finite energy scales. We shall 
further comment on this result below. For frequencies ex- 
ceeding a crossover scale r'o/(27ri?)^, where R is the radius 
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of the nanotube, Egger and Gogolin found a crossover 
to two-dimensionality, which results in a power-law depen- 
dence of the DOS. However, below this crossover scale. Id 
behavior is expected. 

The rest of this work is organized as follows: In Sec. | 
we critically examine the non-perturbative method lead- 
ing to the above results for the zero bias anomaly. In 
particular, we argue that only in the case of long-range 
Coulomb interactions in 2d the method can be formally 
justified; in particular, in quasi Id the approximations 
leading to the result given in Eq. ( p^ ) are uncontrolled, 
so that we have serious doubts whether the exponential 
suppression of the DOS given in Eq. ([l^ ) is the physi- 
cally correct low-frequency behavior in quasi Id. In Sec. 
§ we generalize the calculation of Ref. [0 to finite tem- 
peratures and present numerical results for the frequency- 
dependence of the average DOS which in principle can be 
compared with experimentally measured tunneling con- 
ductances. In Sec. H we discuss the DOS in quasi Id and 
explain why in this case the calculation cannot be per- 
formed with the same confidence as in 2d. Finally, we end 
with a brief summary. 



2 Summing the leading singularities via a 
gauge transformation 

The average DOS of an interacting Fermi system at fi- 
nite temperature T = 1/(3 can be written in terms of 
the disorder-averaged Green function at coinciding space 
points G{i^) = ^(i",r,t^) as 

i^{uj,T) = -^coth (^^^ ImGiuj) . (13) 

To make contact with the disorder-averaged imaginary- 
time Green function at coinciding space points G{t), we 
notice that due to particle-hole symmetry near the Fermi 
energy we have 



Gico) 



nOO 

-2 / dt ainujt G{t 
Jo 



it + 0+) 



(14) 



The problem of calculating the average DOS is now re- 
duced to the problem of calculating the disorder-averaged 
imaginary-time Green function. An attempt to calculate 
within perturbation theory fails at low frequencies 
because the perturbative expansion of v{uj) is plagued by 
singularities, see Eqs. (|l|) and (|^). In 2d with long-range 
Coulomb interactions, these singularities diverge as In^ lo. 
At the same time, however, the perturbative expansion 
of the conductivity cr(w) contains only less severe Incj- 
singularities As a consequence, for sufficiently small 

frequencies, we may re-sum the most singular terms in the 
expansion of the average DOS without considering simul- 
taneously a similar re-summation for the average conduc- 
tivity, because we know a priori that In^ -singularities do 
not appear in the calculation of 17(0;). In this sense, the 
problem of calculating the average DOS decouples from 
the problem of calculating the average conductivity. As 



shown in Ref. [ pO| , under these conditions we may sum 
the most singular contributions to the average DOS via a 
simple gauge transformation, which at finite temperature 
T leads to the following expression for the imaginary time 
Green function, 



G(r) « Go(r)e'3(^) , 
where Go(t) is the Green function of free fermions. 



sin(7rT//3) ' 

and the Debye- Waller factor Q{t) is given by [fsl 



(15) 



(16) 



1 

(3V 



f, 



RPA 



{Doce 



cos w„,,r 



(17) 

Here, f^^^^ is the dynamically screened average Coulomb 
interaction, which is a function of momentum q and 
bosonic Matsubara frequency LUm — 2'inn/ (3. The volume 
of the system is denoted by V . It is important to emphasize 
that in deriving Eqs. ( p"5[]l7| ) the sub-leading logarithmic 
corrections have been ignored, so that it would be incon- 
sistent to retain sub-leading terms involving only a single 
logarithm in the evaluation of the Debye- Waller factor. 

Apparently, the above gauge-transformation method 
to re-sum the leading terms in the perturbative expansion 
of the single-particle Green function has first been used by 
Nazarov (l5|. Later, several authors employed this tech- 
nique to calculate the DOS of disordered interacting elec- 
trons |,|lO[|l|. However, one should keep in mind that in 
practice this method relies on the fact that the problem 
of calculating the average single-particle Green function 
decouples in the sense discussed above from the problem 
of calculating the conductivity. In particular, for 2d dis- 
ordered electrons subject to short-range interactions, the 
perturbative calculation of v{u)) and ct{lij) both involve 
In w-singularities , so that a naive application of the above 
gauge-transformation trick for short-range interactions in 
2d is at least problematic. This is also the case for dif- 
fusive quasi Id electrons, where the perturbative calcula- 
tion of both v{uj) and a{uj) leads to [(^["""^/^-singularities 
[0. Hence, also in this case it is problematic to calculate 
the average DOS using Eqs. (p^|-p7|) without considering 
simultaneously the low-frequency behavior of the conduc- 
tivity. We shall come back to this point in Sec. |4| 



3 Zero bias anomaly in 2d 

In 2d the RPA-interaction is given by [Q 



/n«(2^o^o)-^kn|/q' 



(18) 



This expression is valid in the frequency-momentum 
regime 

\u,n\/Don < |q| « [\lo^\/DoY'^ , M < , (19) 
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which is responsible for the In -corrections to the DOS. 
Converting the sum over Matsubara frequencies into an 
integral, the Debye- Waller factor (analytically continued 
to real time) may be written as 



g(^^) 



■To 



UJ 



— m 

1 — cos{ujt) 
tanh(/3tj/2) 



i sm{u!t) 



(20) 



where we have taken the thermodynamic limit ^ oo. 
To be consistent with the approximation made in Eq. ( p^ , 
we retain only In^-singularities and ignore all terms involv- 
ing only single logarithms: Within this approximation the 
imaginary part of Q{it) vanishes, and for arbitrary tem- 
peratures we find Q{it) ^ ln(t/Ti) ln(t/To). The DOS 

at finite temperature T — 1/P is now given by 



vij-iJ, T) w vq coth 



dt 



sm{ujt) 
sinh(7rt//3) 



X exp 



(21) 



In the limit T — + this equation reduces to Eq. (g) . On the 
other hand, if we first take the limit w — > and then con- 
sider the leading behavior at low temperatures, we obtain 
for T < /22 

K0,T)^2C^, (22) 



where the numerical constant C is given in Eq. (|ll|). Note 
that the prefactor in Eq. (E2) is twice as large as the corre- 
sponding prefactor in Eq. ([lO|). Thus, at low temperatures 
the average DOS at the Fermi energy vanishes linearly in 
T, which resembles the behavior of the average DOS in 
the localized regime. A plot of the DOS for various tem- 
peratures is shown in Fig. ^. 

Formally, the derivation of Eq. ( ^l| ) is only valid in 
the limit of weak disorder where vq — l/irkpi <^ 1. 
However, if the conductivity has a finite DC limit (7(0), 
then it is reasonable to expect that the qualitative be- 
havior of the DOS can be obtained by simply replacing 
r-o ^ n = {e^/h)/TTa{0) in Eq. (|l]). In Rcf. [|o| this re- 
placement has been justified via a simple rcnormalization 
group argument. For a possible comparison of our results 
with future tunneling experiments on semiconductor ma- 
terials which apparently show a metal-insulator transition 
1^,^ (as far as we know, such experiments have not been 
performed yet), we have plotted in Fig. |l| the frequency- 
dependence predicted by Eq. ( ^l|) for tq = I/tt, 
corresponding to the zero-temperature conductivity of or- 
der a{0) « e'^/h. Saturation values in this regime are 
typically encountered in the metallic regime close to the 
apparent metal- insulator transition |^ . The temperature- 
dependence of the average DOS at the Fermi energy is 
shown in Fig. |^. 

It should be noted that in deriving Eq. ( pT| ) we have 
consistently neglected corrections involving single loga- 
rithms. We cannot exclude the possibility that these cor- 
rections lead to a further depletion of the Coulomb gap. 



For example, taking the (sub-leading) imaginary part of 
the full Debye- Waller factor given in Eq. ( pO[ ) into account 
would lead to a nonlinear suppression of the DOS. How- 
ever, as discussed above, such an approximation would not 
be systematic. 




Fig. 1. Graph of the average DOS in 2d for various 

temperatures, see Eq. ([ill). We have chosen ro r, — I/tt 
and Ti — tq. Curves from top to bottom are for T/S72 = 
0.8,0.4,0.2,0.1,0.05,0.025,0. 




Fig. 2. Temperature-dependence of the average DOS i^(T) at 
the Fermi energy in 2d. The parameters are again ro r, = 
I/tt and n = ro . 



4 Zero bias anomaly in Id 

As discussed in Sec. ^, the naive application of Eqs. ( p^ 
in quasi Id is problematic, because in this case the 
interaction corrections to the DOS and to the conductiv- 
ity both involve the same type of singularities. Keeping 
this caveat in mind, let us nevertheless briefly discuss the 



Lorenz Bartosch and Peter Kopietz: Zero bias anomaly in the density of states of low-dimensional metals 



5 





Fig. 3. Graph of the average DOS vit^i) in Id for various tem- 
peratures, see Eq. ( p^ . The curves from top to bottom are for 
T/Qi = 32, 16, 8, 4, 2, 1, 0.5, 0.25, 0. 



Fig. 4. Solid line: graph of the average DOS i^{T) for 
Lu = 0, see Eq. (^^. Dotted dotted line: the approximation 
exp(-v'27Vr). 



predictions of Eqs. (15 -17|) in quasi Id. A similar calcula- 
tion has recently been performed in Ref. [ p^ . 

Since screening is much less effective in one dimen- 
sion than in higher dimensions, we simply approximate the 
one-dimensional RP A- interaction by a constant, f^^^ ~ 
fo 0. This approximation should be correct up to loga- 
rithmic corrections in frequency. A calculation analogous 
to the one leading to Eq. (|20|) results in 



dco 



1 — COs((jjt) 



+ 1 



i\ft 



10 \/2^ tanh(/3u;/2) 

_ (23) 
where fli is given in Eq. (||). Note that in Id, there is 
no need for an ultraviolet cutoff. At T = 0, we have 
Re(5(it) = Im(5(it), such that the imaginary part of 
Q{it) cannot be neglected. The DOS at finite tempera- 
ture T = 1/P can now be written as 



For |a;| <C this integral is easily evaluated within a 
saddle-point approximation, 



i/(w,r = 0) 




exp 



f2i 



\uj\ < I2i , (26) 



which differs from the corresponding expression given by 
Mishchenko et al. (see Eq. (p^)) by a different pref- 
actor. The exponential suppression of the DOS is quite 
surprising and is possibly an artifact of the inconsistent 
exponentiation of the |cj|^^/^-singularities inherent in Eqs. 
@) and (|l|) for quasi Id systems. Recall that for long- 
range Coulomb interactions in 2d it is possible to separate 
In^-terms from sub-leading terms already neglected in the 
approximation made in Eq. (^). In Id, a similar separa- 
tion is not possible. The ultimate low-frequency behavior 
of the DOS in Id might therefore be altered by diagrams 
neglected in Eq. (|l5|). 



vito, T) « vq coth 



X exp 



°° sin(wi) cos(V2]7it) 
sinh(7ri//3) 



duj'- 



1 — cos{Lo't) 
w')3/2 tanh(/3w'/2) 



(24) 



As can be easily checked, for |a;| ^ max{]7i,T|, Eq. ( p4[ ) 
reduces to the perturbative result given in Eq. (0) . A graph 
of v{u>,T) for different temperatures is shown in Fig. ^. 
The DOS at the Fermi energy as a function of temper- 
ature is shown in Fig. ^. The DOS is well approximated 
by v{T) K, vq exp(— ■y/2J7i/r) which is quite similar to a 
result found in Ref. At zero temperature, Eq. ( p4| ) 
may be simplified to 



v{uj,T = 0) = —Re 

TT 



sm X 
dx exp 




(25) 



5 Summary and conclusions 

In summary, we have considered the fate of the singular 
perturbative corrections to the average DOS in quasi Id 
and 2d disordered metals when the frequency and the tem- 
perature are reduced such that one leaves the perturbative 
regime. For a 2d metal with long-range Coulomb interac- 
tions, it is possible to re-sum the leading In^-singularities 
in the perturbative expansion of the average DOS to all 
orders in perturbation theory. The method relies on the 
fact that similar singularities do not appear in the per- 
turbative expansion of the average conductivity. For a 2d 
system with short-range interactions or for a quasi Id sys- 
tem a similar separation of energy scales does not exist, so 
that approximate expressions of the type given in Eq. (p^, 
which involve the exponentiation of a certain subclass of 
Feynman diagrams, are problematic in this case. Thus, in 
practice the gauge-transformation trick described above, 
which has recently been employed by many authors j^. 
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|lO|,|ll, 15|, is only controlled in the case of 2d disordered 
electrons interacting with long-range Coulomb forces in 
the regime where the conductivity has a finite DC-limit. 
In this case this method yields a simple interpolation for- 
mula ( pT[ ) for the average DOS, which predicts a smooth 
crossover from the perturbative regime at high frequen- 
cies to a new low-frequency regime, where u{lo,T) van- 
ishes linearly in w or T. The average DOS of localized 
classical electrons is known to show a similar frequency- 
or temperature-dependence [ p^ . In contrast, the metal- 
lic Coulomb gap discussed here has a quantum mechani- 
cal origin and requires delocalized electrons with a finite 
DC conductivity. Numerical evidence for such a metallic 
Coulomb gap has been found in Ref. . 

In principle, it should be possible to verify the exis- 
tence of the metallic Coulomb gap experimentally via tun- 
neling experiments in strongly correlated disordered sys- 
tems with a finite DC conductivity. The expected shape of 
a typical trace of the tunneling conductance as a function 
of the applied voltage is shown in Fig. |[ Recent tunnel- 
ing experiments by Bielejec et al. ^ in quench-condensed 
Beryllium films show a crossover from the perturbative 
regime with logarithmic corrections to an apparently lin- 
ear Coulomb gap in the DOS. However, at the lowest tem- 
peratures the DOS exhibits a hard correlation gap, the 
origin of which remains open. 
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Actually, in Ref. [ho| the factor Doq^ in the denominator of 
Eq. (|l^) does not explicitly appear; instead, the frequency- 
summation is restricted to the regime luml/Don <g |q| <C 
{\ujm\/ Doy^^ ■ Both prescriptions lead to identical results 
as far as the summation of the leading In'^-singularities 
are concerned. Here we implement the cutoff as in Eq. 
(p^), because then Eqs. ( p^ ) and (jl^) are identical with 
the result obtained by Kamenev and Andreev within a 
field theoretic approach based on the Keldysh formalism 
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